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A new approximate scheme, DSUBm, is described for the coupled cluster method. We then apply it to 
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1. Introduction 

The coupled cluster method (CCM) is a universal microscopic technique of quantum many- 
body theory [ [U [21 El SI [5l [61 [Tl El [3] . It has been applied successfully to many physical systems 
including: 

• systems existing in the spatial continuum, e.g., the electron gas [QUE], atomic nuclei and 
nuclear matter [[El [13], and molecules [Q3]; and 

• systems on a discrete spatial lattice, e.g., spin-lattice models of quantum magnetism [ 1X5 ] fT6 | 

era on lei na nu na lis 12a [iraim. 

A special characteristic of the CCM is that it deals with infinite systems from the outset, and hence 
one never needs to take the limit N — > oo explicitly in the number, N, of interacting particles or the 
number of sites in the lattice. However, approximations in the inherent cluster expansions for the 
correlation operators are required. Several efficient and systematic approximation schemes for the 
CCM have been specifically developed by us previously for use with lattice systems [ [T6ll28ll29ll3"0] . 
Up till now the most favoured and most successful CCM approximation schemes for lattices have 
been the so-called LSUBm and SUBn-m schemes that we describe more fully below in section 03 
Although the LSUBm scheme, in particular, has been shown to be highly successful in practice 
for a wide variety of both frustrated and unfrustrated spin-lattice systems, a disadvantage of this 
scheme is that the number of spin configurations generally rises very rapidly with the truncation 
index m. This motivates us to develop alternative schemes which satisfy one or both of the following 
two criteria: 

• that we are able to calculate more levels of approximation within the scheme, and hence have 
available more data points for the necessary extrapolations for calculated physical quantities 
to the exact limit where all configurations are retained; and 

• that one can capture the physically most important configurations at relatively low orders, 
so that the quantities of interest converge more rapidly with the truncation index. 
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A main aim of our work is thus to provide users of the CCM with more choices of approximation 
scheme. In this paper we outline the formal aspects, of and present preliminary results for some 
benchmark models for, a new CCM approximation scheme, called the DSUBm scheme, for use with 
systems described on a spatial lattice. The rest of the paper is organized as follows. In section [2] 
we introduce the CCM formalism in general. In section [3] we then discuss the specific application 
of the CCM to spin-lattice systems. We consider in section [4] both the existing truncation schemes 
and introduce the alternative new DSUBm scheme. In order to evaluate the accuracy of the new 
approximation scheme we apply it to two very well-studied antiferromagnetic spin-lattice models [ 
[TBI [20], namely: the XX Z and the XY models on the square lattice in two dimensions (2D). Both 
models have a quantum phase transition. They have also been very successfully investigated by 
the CCM within the LSUBm scheme for the ground-state (gs) energy and the gs order parameter 
which is, in our case, the sublattice magnetization. All techniques applied to lattice spin systems 
need to be extrapolated in terms of some appropriate parameters. For exact diagonalization and 
quantum Monte Carlo methods, this is the lattice size N. As noted above, one good aspect of 
the CCM is that we may work in the limit of infinite lattice size (N — > 00) from the outset. By 
contrast, the extrapolation for the CCM is in terms of some truncation index m, where in the 
limit m — > 00 we retain all configurations. The CCM extrapolations that have been used up till 
now in the trunction index, e.g., m for the LSUBm approximation, are heuristic schemes, but we 
have considerable prior experience [ [20] [29] [311 [32l [33] in using and refining them, as described 
in section [5] The new DSUBm approximation scheme is then applied to the square-lattice spin- 
1/2 antiferromagnetic XX Z model in section [6] and the corresondingly XY model in section [7] 
respectively. The results for both models are compared critically with those from the corresponding 
LSUBm scheme and also with the results from other methods. Finally, our conclusions are given 
in section [8] where we reiterate a brief summary of the results. 



2. The CCM Formalism 

This section briefly describes the CCM formalism (and see e.g., Refs. [EJ[9] for further details). A 
first step in any CCM application is to choose a normzalized model (or reference) state |$) that can 
act as a cyclic vector with respect to a complete set of mutually commuting multiconfigurational 
creation operators Cf = {Cjy. The index I here is a set-index that labels the many-particle 
configuration created in the state C7/|$). The exact ket and bra gs energy eigenstates \^f) and 
of the many-body system are then parametrised in the CCM form as: 

Ket-state Bra-state 

|*)=e s |$); (#| = (^e" 5 , (1) 
S = ^<S Z C+; 5=1 + X\£/C7, (2) 

ff|*> = E\V) ; (*\H = E(*\, (3) 

where we have defined Cq = 1 = Cq . The requirements on the multiconfigurational creation oper- 
ators are that any many-particle state can be written exactly and uniquely as a linear combination 
of the states {C/|$)}, which hence fulfill the completeness relation 

£c7|$)($|cy = 1 = |*)(*| + J2ctm(<s>\C7 , (4) 

together with the conditions, 

cy\$) = o = ($\c+; VTyo, (5) 

[C+,C+]=0 = [Cf,Cj]. (6) 
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Approximations are necessary in practice to restrict the label set I to some finite (e.g., LSUBm) 
or infinite (e.g., SUBrt) subset, as described more fully below. The correlation operator S is a linked- 
cluster operator and is decomposed in terms of a complete set of creation operators Cf . It creates 
excitations on the model state by acting on it to produce correlated cluster states. Although the 
manifest Hermiticity, ((*|) t = |*)/(*|#), is lost, the normalization conditions (*|*) = ($|*) = 
($|$) = 1 are preserved. The CCM Schrodinger equations are thus writtern as 

He s \^>) = Ee s \<I>) ; (^\Se~ s H = E(<f>\Se~ s ; (7) 

and its equivalent similarity-transformed form becomes 

e" s He s |$) = E\$) ; (<S>\Se~ s He s = E(<f>\S . (8) 

We note that while the parametrizations of equations ([I]) and J2]) are not manifestly Hermitian 
conjugate, they do preserve the important Hellmann-Feynman theorem at all levels of approxima- 
tions (viz., when the complete set of many-particle configurations I is truncated) [[§]. Furthermore, 
the amplitudes (Si, Si) form canonically conjugate pairs in a time-dependent version of the CCM, 
by contrast with the pairs (Si,Sj) coming from a manifestly Hermitian-conjugate representation 
for = (($|e s e s |$)) _1 ( ( i>|e" s , which are not canonically conjugate to one another. 

The static gs CCM correlation operators S and S contain the real c-number correlation coef- 
ficients Si and Si that need to be calculated. Clearly, once the coefficients {Si, Si} are known, 
all other gs properties of the many-body system can be derived from them. For example, the gs 
expectation value of an arbitrary operator A can be expressed as 

A = (A) = = ($|5e^e s |$) = A(Si,Si) . (9) 

To find the gs correlation coefficients {Si, Si} we simply insert the parametrization of equation 
|(2|) into the Schrodinger equations iJH), and project onto the complete sets of states {{^IC^} and 
{Cf |$)}, respectively, 

($|C7e- 5 #e s |$) = 0; (10) 
{§\S{e- s He s - E)C+\<$>) =0; V/ ^ . (11) 

Equation (fTT]l may also easily be rewritten, by pre-multiplying the ket-state equation |(8|) with the 
state ($>\SCf and using the commutation relation ([6]), in the form 

(<f>\Se- s [H,C+]e s \<S>) = 0; V/ ^ . (12) 

Equations (fT0j) - (fT2"ll may be equivalently derived by requiring that the gs energy expectation value, 
H = (^f\H\^f) = (&\Ser s He s \&), is minimized with respect to the entire set {Si, Si}. In practice 
we thus need to solve equations l|10p and (fT2"ll for the set {Si, Si}. We note that equations |(9]) and 
flOj) show that the gs energy at the stationary point has the simple form 

E = E(Si) = ($\e- s He s \<i>} , (13) 

which also follows immediately from the ket-state equation iJH) by projecting it onto the state (<J>|. 
It is important to note that this (bi-) variational formulation does not necessarily lead to an upper 
bound for E when the summations over the index set {7} for S and S in equation |(2]) are truncated, 
due to the lack of manifest Hermiticity when such approximations are made. Nevertheless, as we 
have pointed out above, one can prove [ |9] that the important Hellmann-Feynman theorem is 
preserved in all such approximations. 

We note that equation (flOf represents a coupled set of multinomial equations for the c-number 
correlation coefficients {Si}. The nested commutator expansion of the similarity-transformed Hamil- 
tonian, 

e- s He s = H+[H,S} + i[[J7, S], S] + ■ ■ ■ , (14) 
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and the fact that all of the individual components of S in the decomposition of equation ^ 
commute with one another by construction [and see equation J6])], together imply that each element 
of S in equation (f2]) is linked directly to the Hamiltonian in each of the terms in equation (f!4|) . Thus, 
each of the coupled equations (fTOf is of Goldstone linked- cluster type. In turn this guarantees that 
all extensive variables, such as the energy, scale linearly with particle number N. Thus, at any level 
of approximation obtained by truncation in the summations on the index / in the parametrizations 
of equation ([2]), we may (and, in practice, do) work from the outset in the limit N — > oo of an 
infinite system. 

Furthermore, each of the seemingly infinite-order (in S) linked-cluster equations lfT0|) will ac- 
tually be of finite length when expanded using equation (fl4|) . since the otherwise infinite series 
of equation lfl4|) will actually terminate at a finite order, provided only (as is usually the case, 
including those for the Hamiltonians considered in this paper) that each term in the Hamiltonian 
H contains a finite number of single-particle destruction operators defined with respect to the 
reference (or generalized vacuum) state |$). Hence, the CCM parametrization naturally leads to 
a workable scheme that can be implemented computationally in an efficient manner to evaluate 
the set of configuration coefficients {Si, Si} by solving the coupled sets of equations (fT0|) and 
lfl2"|). once we have devised practical and systematic truncation hierarchies for limiting the set of 
multiconfigurational set-indices {/} to some suitable finite or infinite subset. 

3. CCM for Spin-Lattice Systems 

We will discuss various practical CCM truncation schemes that fulfill the criteria of being 
systematically improvable in some suitable truncation index to, and that can be extrapolated 
accurately in practice to the exact, m — ► oo, limit for calculated physical quantities. Before doing 
so, however, we first describe how the general CCM formalism described in section[2]is implemented 
for spin-lattice problems in practice. As is the case for any application of the CCM to a general 
quantum many-body system, the first step is to choose a suitable reference state |$) in which the 
the state of the spin (viz., in practice, its projection onto a specific quantization axis in spin space) 
on every lattice site k is characterized. Clearly, the choice of |$) will depend on both the system 
being studied and, more importantly, which of its possible phases is being considered. We describe 
examples of such choices later for the particular models that we utilize here as test cases for our 
new truncation scheme. 

Whatever the choice of |$) we note firstly that is very convenient, in order to set up as universal 
a methodology as possible, to treat the spins on every lattice site in an arbitrarily given model 
state |$) as being equivalent. A suitably simple way of so doing is to introduce a different local 
quantization axis and a correspondingly different set of spin coordinates on each lattice site k, so 
that all spins, whatever their original orientation in |$) in a global spin-coordinate system, align 
along the same direction (which, to be specific, we henceforth choose as the negative z direction) 
in these local spin-coordinate frames. This can always be done in practice by defining a suitable 
rotation in spin space of the global spin coordinates at each lattice site fc. Such rotations are 
canonical transformations that leave unchanged the fundamental spin commutation relations, 

[s+, v] = 2a' k S kk > ; [s z k , 4] = ±s£6 kk , , (15) 

4 = 4±u%, (16) 

among the usual SU(2) spin operators (sjjj, s v k , s|) on lattice site fc. Each spin has a total spin quan- 
tum number, s k , where = Sk{sk + 1) is the SU(2) Casimir operator. For the models considered 
here, s k = s = 1/2, Vfc. 

In the local spin frames defined above the configuration indices I simply become a set of lattice 
site indices, I — > (fci, fc 2 , • • • , fc m ). The corresponding generalized multiconfigurational creation 
operators Cf thus become simple products of single spin-raising operators, Cf — > sjj" sjj" ■ • • s£ , 
and, for example, the ket-state CCM correlation operator is expressed as 

5 = E E s ^ 2 - km si 1 si---st m , a?) 

m k 1 k 2 ---k m 
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and S is similarly defined in terms of the spin-lowering operators sjjT ■ Since the operator S acts on 
the state |$), in which all spins point along the negative z-axis in the local spin-coordinate frames, 
every lattice site fcj in equation (fl7j) can be repeated up to no more than 2s times in each term 
where it is allowed, since a spin s has only (2s + 1) possible projections along the quantization axis. 

In practice the allowed configurations are often further constrained by symmetries in the prob- 
lem and by conservation laws. An example of the latter is provided by the XX Z model considered 
below in section[6l for which we can easily show that the total z-component of spin, sj = Ylk=i s fc' 
in the original global spin coordinates, is a good quantum number since [sj>, H] = in this case. 
Finally, for the quasiclassical magnetically ordered states that we calculate here for both models 
in sections [6] and [3 the order parameter is the sublattice magnetization, M, which is given within 
the local spin coordinates defined above as 

JV N 

M S --W5>^) = --5>|5e-V fe e s |<I>). (18) 

k k 

After the local spin axes have been chosen, the model state thus has all spins pointing downwards 
(i.e., in the negative z-direction, where z is the quantization axis), 

N 

!<£>) = (5?) | l)h ; in the local spin axes, (19) 
fci 

where | J.) = \s, — s) in the usual \s,m s ) notation for single spin states. 

The similarity-transformed Hamiltonian H = e~ s He~ s , and all of the corresponding matrix 
elements in equations ([9|)- (fT3jl and equation lfT8|) . for example, may then be evaluated in the local 
spin coordinate frames by using the nested commutator expansion of equation l[l4)) , the commutator 
relations of equation (fl5l) , and the simple universal relations 

s-|$)=0; Vfc, (20) 

41^)5 Vfc, (21) 
that hold at all lattice sites in the local spin frames. 

4. Approximation schemes 

The CCM formalism is exact if all many-body configurations I are included in the S and S 
operators. In practice, it is necessary to use approximation schemes to truncate the correlation 
operators. 

4.1. Common previous truncation schemes 

The main approximation scheme used to date for continuous systems is the so-called SUBti 
scheme described below. For systems defined on a regular periodic spatial lattice, we have a further 
set of approximation schemes which are based on the discrete nature of the lattice, such as the 
SUBn-m and LSUBn schemes described below. The various schemes and their definitions for 
spin-lattice systems are: 

1. the SUBrt scheme, in which only the correlations involving n or fewer spin-raising operators 
for S are retained, but with no further restrictions on the spatial separations of the spins 
involved in the configurations; 

2. the SUBn-m scheme wich includes only the subset of all n-spin-flip configurations in the 
SUBn scheme that are defined over all lattice animals of size < m, where a lattice animal is 
defined as a set of contiguous lattice sites, each of which is nearest-neighbour to at least one 
other in the set; and 
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3. the LSUBm scheme which includes all possible multi-spin-flip configurations defined over all 
lattice animals of size < m. The LSUBm scheme is thus equivalent to the SUBn-m scheme 
for n = 2sm for particles of spin quantum number s. For example, for spin-1/2 systems, for 
which no more than one spin-raising operator, s^T, can be applied at each site k, LSUBm = 
SUBm-m. 



4.2. The new DSUBm scheme 

Our new DSUBm scheme is now defined to include in the correlation operator S all possible 
configurations of spins involving spin-raising operators where the maximum length or distance of 
any two spins apart is defined by L m , where L m is a vector joining sites on the lattice and the 
index m labels lattice vectors in order of size. Hence DSUB1 includes only nearest-neighbours, etc. 

Table [TJ shows how L m progresses in terms of k and I (which are the sides of the lattice in the x 
and y directions) for the case of a 2D square lattice with sides parallel to the x and y axes. Similar 
tables can be constructed for an arbitrary regular lattice in any number of dimensions. Table [TJ 
shows, for example, that the DSUB5 approximation on a 2D square lattice involves all clusters of 
spins (and their associated spin-raising operators) for which the maximum distance between any 
two spins is v8- Clearly the DSUBm scheme orders the multispin configurations in terms, roughly, 
of their compactness, whereas the LSUBm scheme orders them according to the overall size of the 
lattice animals (or polyominos), defined as the number of contiguous lattice sites involved. 



5. Extrapolation schemes 

Any of the above truncated approximations clearly becomes exact when all possible multispin 
cluster configurations are retained, i.e., in the limit as n — ► oo and/or m — * oo. We have considerable 
experience, for example, with the appropriate extrapolations for the LSUBm scheme [ 20, 29, 31, 32, 
[33] , that shows that the gs energy behaves in the large-m limit as a power series in 1/m 2 , whereas 
the order parameter M behaves as a power series in 1/m (at least for relatively unfrustrated 
systems). Initial experience with the new DSUBm scheme shows, perhaps not unsurprisingly, that 
in the corresponding large m limit the gs energy and order parameter behave as power series in 
1/i^j and l/L m , respectively, as we show in more detail (below) for the two examples of the spin- 
1/2 XX Z and XY models on the 2D square lattice. It is clear on physical grounds that the index 
L m should provide a better extrapolation variable than the index m itself for the DSUBm scheme, 
and so it turns out in practice. For the present, where we are interested primarily in a preliminary 
investigation of the power and accuracy of the DSUBm scheme, we limit ourselves to retaining 
only the leading terms in the power-series expansions, 



E 
N 



DSUBm 



°o + oi ( TT ) > ( 22 ) 



Table 1. The formulation of the length parameter L m of the DSUBm approximation on a square 
lattice, in terms of lattice increments k and I along the two sides of the lattice. 
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M 



DSUBm 




(23) 



Further sub-leading terms in each of the power series can easily be retained later should it prove 
useful to increase the accuracy of the extrapolations. 

5.1. Three fundamental rules for the selection and extrapolation of the CCM raw data 

We list below three fundamental rules as guidelines for the selection and extrapolation of the 
CCM raw data, using any approximation scheme: 

• RULE 1: In order to fit well to any fitting formula that contains n unknown parameters, one 
should have at least (n + 1) data points. This rule takes precedence over all other rules, and 
is vital to obtain a robust and stable fit. 

• RULE 2: Avoid using the lowest data points (e.g., LSUB2, SUB2-2, DSUB1, etc.) wherever 
possible, since these points are rather far from the large-m limit, unless it is necessary to do 
so to avoid breaking RULE 1, e.g., when only n data points are available. 

• RULE 3: If RULE 2 has been broken (e.g., by including LSUB2 or SUB2-2 data points), 
then do some other careful consistency checks on the robustness and accuracy of the results. 

6. The spin-1/2 antiferromagnetic XXZ model on the square lattice 

In this section, we shall consider the spin-1/2 XXZ model on the infinite square lattice. The 
Hamiltonian of the XXZ model, in global spin coordinates, is written as 



where the sum on (i, j) runs over all nearest-neighbour pairs of sites on the lattice and counts 
each pair only once. Since the square lattice is bipartite, we consider N to be even, so that each 
sublattice contains (1/2)N spins, and we consider only the case where N — > oo. The Neel state is 
the ground state (GS) in the trivial Ising limit A — * oo, and a phase transition occurs at (or near 
to) A = 1. Indeed, the classical GS demonstrates perfect Neel order in the z-direction for A > 1, 
and a similar perfectly ordered x-y planar Neel phase for —1 < A < 1. For A < — 1 the classical 
GS is a ferromagnet. 

The case A = 1 is equivalent to the isotropic Heisenberg model, whereas A = is equivalent 
to the isotropic version of the XY model considered in section [7| below. The z component of total 
spin, sij,, is a good quantum number as it commutes with the Hamiltonian of equation l|24p . Thus 
one may readily check that [sf>, Hxxz] = 0. Our interest here is in those values of A for which the 
GS is an antiferromagnet. 

The CCM treatment of any spin system is started by choosing an appropriate model state |<fr) 
(for a particular regime) , so that a linear combinations of products of spin-raising operators can be 
applied to this state and all possible spin configurations are determined. There is never a unique 
choice of model state |<f>). Our choice should clearly be guided by any physical insight that we 
can bring to bear on the system or, more specifically, to that particular phase of it that is under 
consideration. In the absence of any other insight into the quantum many-body system it is common 
to be guided by the behaviour of the corresponding classical system (i.e., equivalently, the system 
when the spin quantum number s — > oo). The XXZ model under consideration provides just such 
an illustrative example. Thus, for A > 1 the classical Hamiltonian of equation (|24| on the 2D square 
lattice (and, indeed, on any bipartite lattice) is minimized by a perfectly antiferromagnetically 
Neel-ordered state in the spin z-direction. However, the classical gs energy is minimized by a Neel- 
ordered state with spins pointing along any direction in the spin x-y plane (say, along the spin 
x-direction) for — 1 < A < 1. Either of these states could be used as a CCM model state |<£>) 
and both are likely to be of value in different regimes of A appropriate to the particular quantum 




(24) 



(id) 
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phases that mimic the corresponding classical phases. For present illustrative purposes we restrict 
ourselves to the z-aligned Neel state as our choice for |$), written schematically as 

|<3>) = | ■ ■ • HIT ' ' ') i m the global spin axes, 

where | |) = ||, +5) and | J.) = |, — 1) in the usual \s,m s ) notation. Such a state is, clearly, likely 
to be a good starting-point for all A > 1, down to the expected phase transition at A = 1 from a 
z-aligned Neel phase to an x-y planar Neel phase. 

As indicated in section [3] it is now convenient to perform a rotation of the axes for the up- 
pointing spins (i.e., those on the sublattice with spins in the positive z-direction) by 180° about 
the spin y-axis, so that |<I>) takes the form given by equation lfl9|) . Under this rotation, the spin 
operators on the original up sub-lattice are transformed as 

s x -s x , s v s y , s z ™> ~s z . (25) 

The Hamiltonian of equation (|24| may thus be rewritten in these local spin coordinate axes as 

Hxxz = -\ i44 + S ^ S J + 2As *4 ■ (26) 

(id) 

As in any application of the CCM to spin-lattice systems, we now include in our approxima- 
tions at any given order only those fundamental configurations that are distinct under the point 
and space group symmetries of both the lattice and the Hamiltonian. The number, Nf, of such 
fundamental configurations at any level of approximation may be further restricted whenever ad- 
ditional conservation laws come into play. For example, in our present case, the XX Z Hamiltonian 
of equation (|24| commutes with the total uniform magnetization, — X)fe=i s fc; m the global spin 
coordinates, where k runs over all lattice sites. The GS is known to lie in the — subspace, and 
hence we exclude configurations with an odd number of spins or with unequal numbers of spins on 
the two equivalent sublattices of the bipartite square lattice. We show in figure Q] the fundamental 
configurations that are accordingly allowed for the DSUBm approximations for this spin-1/2 XX Z 
model on the 2D square lattice, with 1 < m < 4. We see that, for example, Nf — 12 at the DSUB4 
level of approximation. We also see that Nf = 2 for both the DSUB2 and DSUB3 approximations, 
since the conservation law = does not permit any additional configurations of spins with a 
maximum distance L3 — 2, apart from those already included in the DSUB2 approximation. 



6.1. Ground-state energy and sublattice magnetization 

The DSUBm approximations can readily be implemented for the present spin-1/2 XX Z model 
on the 2D square lattice for all values m < 11 with reasonably modest computing power. By 
comparison, the LSUBm scheme can be implemented with comparable computing resources for all 
values m < 10. Numerical results for the gs energy per spin and sublattice magnetization are shown 
in table[2]at the isotropic point A = 1 at various levels of approximation, and corresponding results 
for the same quantities are displayed graphically in figures [2] and [3] as functions of the anisotropy 
parameter A. 

We also show in table [2] for the isotropic Heisenberg Hamiltonian (A = 1) the results for the gs 
energy and sublattice magnetization using the leading (linear) extrapolation schemes of equations 
(|22|) and ((23)1 respectively of the DSUBm data, employing various subsets of results. Comparison 
is also made with corresponding LSUBm extrapolation schemes for the same model [ [3U I35j . 
The results are generally observed to be in excellent agreement with each other, even though 
the DSUBoo extrapolations have employed the simple leading (linear) fits of equations (|22)l and 
(|23|) . whereas the corresponding LSUBoo results shown [ [3U [35] have been obtained from the 
potentially more accurate quadratic fits E/N = a + ai(l/m 2 ) + a 2 (l/m 2 ) 2 , M = b + &i(l/m 2 ) + 
6 2 (l/w) 2 , to the LSUBm data, in which the next-order (quadratic) corrections have also been 
included in the relevant expansion parameters, 1/m and 1/m, respectively. Excellent agreement 
of all the CCM extrapolations is also obtained with the results from the best of the alternative 
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Figure 1. The fundamental configurations for the DSUBm scheme with m = {1,2,3,4} for the 
spin-1/2 XX Z model on a square lattice in two dimensions. The filled circles mark the relative 
positions of the sites of the square lattice on which the spins are flipped with respect to the 
model state. The unfilled circles represent unflipped sites. 



methods for this model, including third-order spin- wave theory (SWT) [[36], linked-cluster series 
expansion techniques [ 38|, and the extrapolations to infinite lattice size (iV — > oo) from the exact 
diagonalization (ED) of small lattices [[37], and quantum Monte Carlo (QMC) calculations for 



DSUBl 

DSUB2=DSUB3 

DSUB4 

DSUB6 

DSUB8 

DSUB10 

DSUBoo 




4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 



Figure 2. (Color online) CCM results for the ground-state energy per spin, E/N, as a function 
of the anisotropy parameter A, of the spin-1/2 XX Z model on the square lattice, using various 
DSUBm approximations based on the ^-aligned Neel model state. The DSUBm results with 
m = {6, 8, 10} are extrapolated using the leading (linear) fit of equation ll22l and shown as the 
curve DSUBoo. Ai = magnetization point of inflexion, described in the text. 
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Table 2. The ground-state energy per spin (E/N) and sublattice magnetization (M) for the 
spin-1/2 XXZ model on the square lattice, obtained using the CCM DSUBm approximation 
scheme with l<m<llatA = l. Nf is the number of fundamental configurations at a given 
DSUBm or LSUBm level of approximation. Ai = DSUBm sublattice magnetization point of 
inflexion. The DSUBm results for odd values of m, even values of m and the whole series of m 
are extrapolated separately. These results are compared to calculations using third-order spin- 
wave theory (SWT), exact diagonalization (ED), series expansion (SE), quantum Monte Carlo 
(QMC) and LSUBoo extrapolations of the CCM LSUBm approximations. 



Method 




N f 


E/N 


M 


A, 


Max. No. 








A = l 




of spins 


DSUB1=LSUB2 


1 


1 


-0.64833 


0.421 





2 


DSUB2 


V2 


2 


-0.65311 


0.410 


0.258 


4 


DSUB3 


2 


2 


-0.65311 


0.410 


0.258 


4 


DSUB4 


V5 


12 


-0.66258 


0.385 


0.392 


6 


DSUB5 


V8 


20 


-0.66307 


0.382 


0.479 


8 


DSUB6 


3 


43 


-0.66511 


0.375 


0.506 


8 


DSUB7 


VTo 


135 


-0.66589 


0.371 


0.629 


12 


DSUB8 




831 


-0.66704 


0.363 


0.614 


14 


DSUB9 


4 


1225 


-0.66701 


0.363 


0.654 


14 


DSUB10 


VU 


6874 


-0.66774 


0.357 


0.637 


16 


DSUB11 


Vl8 


14084 


-0.66785 


0.356 


a 


16 


LSUB8 




1287 


-0.66817 


0.352 




8 


LSUB10 




29605 


-0.66870 


0.345 




10 




Extrapolations 




Based on 


E/N 


M 


A, 


DSUBoo 


m = 


{6,8,10} 


-0.67082 


0.308 


1.009 


DSUBoo 


m = {5,7,9,11} 


-0.67122 


0.311 


1.059 b 


DSUBoo 


m = 


{7,9,11} 


-0.66978 


0.319 




DSUBoo 


3 < 


m < 11 


-0.67177 


0.315 


1.025 c 


DSUBoo 


4 < m < 11 


-0.66967 


0.325 


0.979 


LSUBoo [ 34J 


m — 


{3,5,7,9} 


-0.67029 


0.305 




LSUBoo [SUED 


m= {4,6,8,10} 


-0.66966 


0.310 




LSUBoo 


m = 


{6,8,10} 


-0.66962 


0.308 




SWT [[36] 






-0.66999 


0.307 




SE [ 38J 






-0.6693 


0.307 




ED [[37] 






-0.6700 


0.3173 




QMC [[39] 






-0.669437(5) 


0.3070(3) 





NOTES : 

a The magnetization point of inflexion for DSUB11 is not available since we only calculated at 
A = 1 in this approximation. 

b The magnetization points of inflexion for the odd DSUBm levels are extrapolated using m = 
{5,7,9}. 

c The magnetization points of inflexion for the whole series of DSUBm data are extrapolated as 
indicated, but without m = 11. 

larger lattices [ 40J . 

We note that it has been observed and well documented in the past (and see, e.g., Ref. [ 34J) 
that the CCM LSUBm results for this model (and many others) for both the gs energy E and the 
sublattice magnetization M show a distinct period-2 "staggering" effect with index m, according 
to whether m is even or odd. As a consequency the LSUBm data for both E and M converge 
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Figure 3. (Color online) CCM results for the sublattice magnetization, M, as a function of the 
anisotropy parmeter A, of the spin-1/2 XXZ model on the square lattice, using various DSUBm 
approximations based on the ^-aligned Neel model state. The DSUBm results with m = {6, 8, 10} 
are extrapolated using the leading (linear) fit of equation Il23ll and shown as the curve DSUBoo. 
Aj = point of inflexion in the curve, shown by arrows in the figure. 



differently for the even-m and the odd-m sequences, similar to what is observed very frequently in 
perturbation theory in corresponding even and odd orders [[40]. As a rule, therefore, the LSUBm 
data are generally extrapolated separately for even m and for odd values of m, since the staggering 
makes extrapolations using both odd and even values together extremely difficult. We show in 
figure H] our DSUBm results for the gs energy per spin and the sublattice magnetization plotted 
against and 1/L m , respectively. The higher m values clearly cluster well in both cases on 

straight lines, thereby justifying a posteriori our heuristic extrapolation fits of equations f22|) and 
lj23|) . Just as in the LSUBm case a slight "even-odd staggering" effect is observed in the DSUBm 
data (perhaps more so for the energy than for the sublattice magnetization), although it is less 
pronounced than for the corresponding LSUBm data [134]. 

6.2. Termination or critical points 

Before discussing our DSUBm results further for this model we note that the comparable 
LSUBm solutions actually terminate at a critical value A c = A c (m), which itself depends on the 
truncation index m [[30]. Such LSUBm termination points are very common for many spin-lattice 
systems and are very well documented and understood (and see, e.g., Ref. [ 30]). In all such cases 
a termination point always arises due to the solution of the CCM equations becoming complex 
at this point, beyond which there exist two branches of entirely unphysical complex conjugate 
solutions [[30]. In the region where the solution reflecting the true physical solution is real there 
actually also exists another (unstable) real solution. However, only the (shown) upper branch of 
these two solutions reflects the true (stable) physical GS, whereas the lower branch does not. The 
physical branch is usually easily identified in practice as the one which becomes exact in some 
known (e.g., perturbative) limit. This physical branch then meets the corresponding unphysical 
branch at some termination point beyond which no real solutions exist. The LSUBm termination 
points are themselves also reflections of the quantum phase transitions in the real system, and may 
be used to estimate the position of the phase boundary [130]. 

It is interesting and intriguing to note that when the DSUBm approximations are applied to 
the XXZ model, they do not terminate as do the corresponding LSUBm approximations. We 
have no real explanation for this rather striking difference in behaviour for two apparently similar 
schemes applied to the same model. However, it is still possible to use our DSUBm data to extract 



11 



R.F. Bishop et al. 
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(a) Ground-state energy per spin 
0.55 | , , , , 1 




0.3 ' 1 1 ' 
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(b) Ground-state sublattice magnetization 

Figure 4. (Color online) Illustration of the staggered nature of the DSUBm scheme for the gs 
energy per spin, E/N, and sublattice magnetization, M, for the spin-1/2 antiferromagnetic 
XX Z model on the square lattice. The DSUBm data are plotted against for E/N and 

against 1/L m for M. The results clearly justify the heuristic extrapolation schemes of equations 
lf22j) and lt23ll . 



an estimate for the physical phase transition point at which the ^-aligned Neel phase terminates. 
As has been justified and utilized elsewhere [ 18J, a point of inflexion at A = Aj in the sublattice 
magnetization M as a function of A clearly indicates the onset of an instability in the system. 
Such inflexion points A< = Aj(m) occur for all DSUBm approximations, as indicated in table [2] 
and figure [31 The DSUBm approximations are thus expected to be unphysical for A < Aj(m), 
and we hence show the corresponding results for the gs energy per spin in figure [2] only for values 
A, > Aj(m), Heuristically, we find that the magnetization inflexion points Aj(m) scale linearly 
with 1/L m to leading order, and the extrapolated results shown in table [2] have been performed 
with the leading (linear) fit, Aj = c Q + c±(l/L m ), commensurate with the corresponding linear 
fits in (1/i^j) and (1/L m ) for the gs energy per spin and sublattice magnetization of equations 
|22|) and i[23|) , respectively. All of the various extrapolations shown in table [2] for Aj in the limit 
m — > oo are in good agreement with one another, thereby again demonstrating the robust quality 
of the heuristic extrapolation scheme. Furthermore, they are also in excellent agreement with the 
expected phase transition point at A c = 1 between two quasiclassical Neel-ordered phases aligned 
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along the spin z-axis (for A > 1) and in some arbitrary direction in the spin x-y-plane (for | A| < 1). 

Although we do not do so here, the x-y planar Neel phase could itself also easily be investigated 
by another CCM DSUBm series of calculations based on a model state |<&) with perfect Neel 
ordering in, say, the x-direction. 

Summarizing our results so far, we observe that the DSUBm scheme has at, least partially, 
fulfilled the expectations placed upon it for the present model. Accordingly, we now apply it to the 
second test model of the spin-1/2 XY model on the 2D square lattice. 



7. The spin-1/2 XY model on the square lattice 

The Hamiltonian of the XY model [ 18J in global spin coordinates, is written as 

fl^ = 3(l+AK^ + (l-A)^; -1<A<1, (27) 

where the sum on again runs over all nearest-neighbour pairs of lattice sites and counts each 
pair only once. We again consider the case of spin-1/2 particles on each site of an infinite square 
lattice. 

For the classical model described by equation (|2"T)l . it is clear that the GS is a Neel state in the 
^-direction for < A < 1 and a Neel state in the y-direction for —1 < A < 0. Hence, since we only 
consider the case < A < 1, we choose as our CCM model state |$) for the quantum XY model 
a Neel state aligned along the x-direction, written schematically as, 

|$) = | • • • < >< ►■••), in the global spin axes. 

Clearly the case — 1 < A < is readily obtained from the case < A < 1 by interchange of the x- 
and y-axes. 

Once again we now perform our usual rotation of the spin axes on each lattice site so that |$) 
takes the form given by equation 1(19)1 in the rotated local spin coordinate frame. Thus, for the 
spins on the sublattice where they point in the negative ^-direction in the global spin axes (i.e., the 
left-pointing spins) we perform a rotation of the spin axes by +90° about the spin y-axis. Similarly, 
for the spins on the other sublattice where they point in the positive x-direction in the global spin 
axes (i.e., the right-pointing spins) we perform a rotation of the spin axes by —90° about the spin 
y-axis. Under these rotations the spin operators are transformed as 

s x — > s z , s v — > s v , s z — > — s x , left-pointing spins; (28a) 
s x — ► -s z , s y — ► s y , s z — > s x , right-pointing spins. (28b) 

The Hamiltonian of equation (|27| may thus be rewritten in the local spin coordinate axes 
defined above as 



Hxy = 

(i,3) 



-(1 + A) S f S | - \{l A)(s+s+ + s-sj) + 1(1 - A)(s+ S - 



(29) 



As before, we now have to evaluate the fundamental configurations that are retained in the 
CCM correlation operators S and S at each DSUBm level of approximation. Although the point 
and space group symmeries of the square lattice (common to both the XX Z and XY models 
considered here) and the two Hamiltonians of equations ([261 and (|29j) are identical, the numbers 
Nf of fundamental configurations for a given DSUBm level are now larger (except for the case 
m = 1) for the XY model than for the XX Z model, since the uniform magnetization is no longer 
a good quantum number for the XY model, [Hxy, Sj] ^ 0. Nevertheless, we note from the form of 
equation l|29p. in which the spin-raising and spin- lowering operators appear only in combinations 
that either raise or lower the number of spin flips by two (viz., the sf Sj~ and s^sj combinations, 
respectively) or leave them unchanged (viz., the s^sj and s~s^ combinations), it is only necessary 
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Figure 5. The fundamental configurations for the DSUBm scheme with m = {1,2,3} for the 
spin-1/2 XY model on a square lattice in two dimensions. The filled circles mark the relative 
positions of the sites of the square lattice on which the spins are flipped with respect to the 
model state. The unfilled circles represent unflipped sites. 



for the Sj. = GS to consider fundamental configurations that contain an even number of spins. 
Thus, the main difference for the XY model over the XXZ model is that we must now also 
consider fundamental configurations in which we drop the restriction for the former case of having 
an equal number of spins on the two equivalent sublattices of the bipartite square lattice that was 
appropriate for the latter case. We show in figure[5]the fundamental configurations that are allowed 
for the spin-1/2 XY model on the square lattice for the DSUBm approximation with 1 < m < 3, 
and we invite the reader to compare with the corresponding fundamental configurations for the 
spin-1/2 XXZ model on the same square lattice shown in figure 1. The corresponding numbers Nf 
of fundamental configurations for the XY model are also shown in table [3] for the higher DSUBm 
approximations with m < 9 for which we present results below. 

7.1. Ground-state energy and sublattice magnetization 

We present results for the spin-1/2 XY model on the square lattice in the CCM DSUBm 
approximations for all values m < 9 that can be easily computed with very modest computing 
power. Comparable computing power enables the corresponding LSUBm scheme to be implemented 
for all m < 8. Numerical results for the gs energy per spin and sublattice magnetization are shown 
in table [3] at the isotropic point at A = at various levels of approximation, and corresponding 
results for the same gs quantities are shown graphically in figures [6] and [7] as functions of the 
anisotropy parameter A. 

We also show in table [3] for the isotropic XY Hamiltonian (A = 0) the results for the gs 
energy and sublattice magnetization using the leading (linear) extrapolation schemes of equations 
|22p and (|23|) respectively of the DSUBm data, employing various subsets of our results, as for 
the XXZ model considered previously. We also compare in table [3] the present results with the 
corresponding CCM LSUBm results [[18] for the same model. All of the CCM results are clearly in 
excellent agreement both with one another and with the results of best of the alternative methods 
available for this model, including the linked-cluster series expansion (SE) techniques [[41] and a 
quantum Monte Carlo (QMC) method [142]. 

We again show in figure [8] our DSUBm results for the present XY model for the gs energy per 
spin and the sublattice magnetization, plotted respectively against l/I^i an ^ V-^m- As previously 
for the XXZ model, the higher m values cluster well on straight lines in both cases, thereby 
justifying once more our heursitic choice of extrapolation fits indicated in equations l]22H and (|23|1 . 
Figures 8(a) and |8(b)| again show an "even-odd" staggering effect in the termination index m for the 
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Table 3. The ground-state energy per spin (E/N) and sublattice magnetization (M) for the spin- 
1/2 XY model on the square lattice, obtained using the CCM DSUBm approximation scheme 
with 1 < m < 9 at A = 0. Nf is the number of fundamental configurations at a given level of 
DSUBm or LSUBm approximation. A c = DSUBm termination point. The DSUBm results for 
odd values of m, even values of m and the whole series of m are extrapolated separately. These 
results are compared to calculations of series expansion (SE), quantum Monte Carlo (QMC) 
and LSUBoo extrapolations of the CCM LSUBm approximations. 



Method 




N f 


E/N 


M 


A c 


Max. No. 








A = 




of spins 


DSUB1=LSUB2 


1 


1 


-0.54031 


0.475 




2 


DSUB2 


v/2 


3 


-0.54425 


0.467 




4 


DSUB3 


2 


6 


-0.54544 


0.464 




4 


DSUB4 


v/5 


21 


-0.54724 


0.458 


-0.253 


6 


DSUB5 


V8 


44 


-0.54747 


0.456 


-0.205 


8 


DSUB6 


3 


78 


-0.54774 


0.455 


-0.181 


8 


DSUB7 


yio 


388 


-0.54811 


0.453 


-0.135 


12 


DSUB8 


v/13 


1948 


-0.54829 


0.451 


-0.107 


14 


DSUB9 


4 


3315 


-0.54833 


0.451 


-0.099 


14 


LSUB6 




131 


-0.54833 


0.451 


-0.073 


6 


LSUB8 




2793 


-0.54862 


0.447 


-0.04 


8 




Extrapolations 




Based on 


E/N 


M 


A c 


DSUBoo 


m = {4,6,8} 


-0.54879 


0.442 


-0.036 


DSUBoo 


m = {5,7,9} 


-0.54923 


0.437 


0.011 


DSUBoo 


4 < 


m < 9 


-0.54884 


0.441 


-0.029 


LSUBoo [ 18J 


m = {4,6,8} 


-0.54892 


0.435 


0.00 


SE [ 41J 






-0.54883 


0.43548 


0.0 


QMC [B2] 






-0.548824(2) 


0.437(2) 





DSUBm data, which is perhaps slightly more pronounced than that for the XX Z model shown in 
figures 4(a) and |4(b)| For this reason we have again shown separate extrapolations of our DSUBm 
results in table [3] for the even-m data and the odd-m data, as well as results using all (higher) 
values of m. 



7.2. Termination or critical points 

It is interesting to note that for the present XY model the CCM DSUBm solutions (with our 
choice of model state as a Neel state in the ^-direction) now do physically terminate for all values 
of the truncation index m > 4 at a critical value A c = A c (m), exactly as commonly occurs (as for 
the present model) for the LSUBm calculations, as we explained above in section [621 Why such 
DSUBm terminations occur for the XY model but not for the previous XX Z model is not obvious 
to us. The corresponding termination points, A c = A c (m), at various DSUBm and LSUBm levels 
of approximation are shown in table [3l It has been shown previously [ 29J that A c (m) scales well 
with 1/m 2 for the LSUBm data, and the LSUBoo result [[18] shown in table [3] was obtained by a 
leading (linear) fit, A c (m) = d + di{\/m 2 ). We find heuristically that the best large-m asymptotic 
behaviour of the DSUBm data for A c (m) is against as the scaling parameter. Accordingly, the 

DSUBoo values for A c in tableware obtained with the leading (linear) fit, A c (m) = do + di(l/L^ n ). 
We see that both the LSUBoo and DSUBoo results for A c = A c (oo) agree very well with the value 
A c = that is known to be the correct value for the phase transition in the one-dimensional 
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Figure 6. (Color online) CCM results for the gs energy of the spin-1/2 XY model on the square 
lattice obtained using the DSUBm approximation based on the Neel state aligned along any 
axis in the x-y plane. The DSUBm results with m = {5, 7, 9} are extrapolated using equation 
|T2"2"|) to give the curve labelled DSUBoo. 




Figure 7. (Color online) CCM results for the gs sublattice magnetization of the spin-1/2 XY 
model on the square lattice obtained using various DSUBm approximations based on the Neel 
state aligned along any axis in the x-y plane. The DSUBm results with m = {5, 7, 9} are 
extrapolated using equation lf23ll to give the curve labelled DSUBoo. 



spin-1/2 XY chain from the known exact solution [[43], and which is believed also to be the 
phase transition point for higher dimensions, including the present 2D square lattice, on symmetry 
grounds. 

8. Conclusions 

From the two nontrivial benchmark spin-lattice problems that we have investigated here, it is 
clear that the new DSUBm approximation scheme works well for calculating their gs properties 
and phase boundaries. We have utilized here only the simplest leading-order extrapolation schemes 
in the pertinent scaling variables, and have shown that these may be chosen, for example, as 1/i^ 
for the gs energy and l/L m for the order parameter. Clearly, in general, the results can be further 
improved by keeping higher-order terms in these asymptotic expansions (i.e., by retaining higher 
powers in the polynomial scaling expansions) although more data points may then be needed, 
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Figure 8. (Color online) Illustration of the staggered nature of the DSUBm scheme for the gs 
energy per spin, E/N, and sublattice magnetization, M, for the spin-1/2 XY model on the 
square lattice. The DSUBm data are plotted against for E/N and against 1/L m for M. 

The results clearly justify the heuristic extrapolation schemes of equations (1221) and (12311 . 



especially in cases where the "even-odd" staggering effect is pronounced, as for XY model presented 
here. For further use of the scheme for more complex lattice models (e.g., those exhibiting geometric 
or dynamic frustration) it will be necessary to re- visit the validity of these expansions, but a great 
deal of previous experience in such cases for the LSUBm scheme will provide good guidance. 

On the basis of the test results presented here, the DSUBm scheme clearly fulfills the first 
of our two main criteria for introducing it, viz., that the number of fundamental configurations, 
Nf, increases less rapidly with truncation index m than for the corresponding LSUBm series of 
approximations. At the same time our second criterion of capturing the physically most important 
configurations at relatively low levels of approximation also seems to be fulfilled, according to our 
experience with the convergence of the DSUBm sequences for observable quantities. At the very 
least we now have two schemes (LSUBm and DSUBm) available to us for future investigations, 
each of which has its own merits, and which thus allows us more freedom in applications of the 
CCM to other spin-lattice models in future. 

The one slight drawback in the scheme which mitigates against our goal of obtaining more 
DSUBm data points, for the same computing power than for the LSUBm scheme applied to the 
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same system, and that hence can be used together to attain more accuracy in the extrapolations, 
is the slight "even-odd" staggering in the data that is observed in the DSUBto results, albeit that 
it is somewhat reduced from the similar stagerring in the corresponding LSUBm results. We have 
some ideas on how the DSUBm scheme might itself be modified to reduce this staggering and we 
hope to report results of these further investigations in a future paper. 
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